Lifetime of S and a limit for its 2p decay 

energy 

A.S. Fomichev/I.G. Mukha,^'^ S.V. Stepantsov/ 
L.V. Grigorenko, ^'^''^ E.V. Litvinova, V. Chudoba, ^'^ 
LA. Egorova,g M.S. Golovkov/ A.V. Gorshkov, ^ 
V.A. Gorshkov, ^ G. Kaminski, ^'^ S.A. Kmpko, ^ 
Yu.L. Parfenova, ^'^ S.I. Sidorchuk, ^ R.S. Slepnev, ^ 
G.M. Ter-Akopian/ R. Wolski/'^ M.V. ZhukovJ 

^Flerov Laboratory of Nuclear Reactions, JINR, RU-141980 Dubna, Russia 

^ GSI Helmholtzzentrum fiir Schwerionenforschung, Planckstrafie 1, D-64-291 

Darmstadt, Germany 

^ CSIC-IFIC, University of Valencia, E-46071 Valencia, Spain 

'^Russian Research Center "The Kurchatov Institute", Kurchatov sq. 1, 
RU- 123 182 Moscow, Russia 

^Institut fur Theoretische Physik, Goethe- Universitdt, 60438 Frankfurt am Main, 

Germany 

^Institute of Physics, Silesian University in Opava, Bezrucovo nam. 13, 74601 

Gzech Republic 

^Bogolubov Laboratory of Theoretical Physics, JINR, Dubna, RU-141980 Russia 

^Institute of Nuclear Physics PAN, Radzikowskiego 152, PL-31342 Krakow, 

Poland 

^Skobeltsyn Institute of Nuclear Physics, Moscow State University, 119991 

Moscow, Russia 

^Fundamental Physics, Ghalmers University of Technology, S-41296 Goteborg, 

Sweden 



Abstract 

Unknown isotope ^^S, expected to decay by two-proton {2p) emission, was studied 
theoretically and was searched experimentally. The structure of this nucleus was 
examined within the relativistic mean field (RMF) approach. A method for taking 
into account the many-body structure in the three-body decay calculations was 
developed. The results of the RMF calculations were used as an input for the three- 
cluster decay model worked out to study a possible 2p decay branch of this nucleus. 
The experimental search for ^^S was performed in fragmentation reactions of a 50.3 
A MeV ^^S beam. No events of ^^S or -^^P (a presumably proton- unstable subsystem 
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of26S) were observed. Based on the obtained production systematics an upper half- 
life limit of Tx/2 < 79 ns was established from the time-of- flight through the fragment 
separator. Together with the theoretical lifetime estimates for two-proton decay this 
gives a decay energy limit of Q2p > 640 keV for ^^S. Analogous limits for ^^P are 
found as T1/2 < 38 ns and Qp > 110 keV. In the case that the one-proton emission 
is the main branch of the ^^S decay a limit Q2p > 230 keV would follow for this 
nucleus. It is likely that ^^S resides in the picosecond lifetime range and the further 
search for this isotope is prospective for the decay-in-fiight technique. 

Key words: two-proton radioactivity, relativistic mean field, three-cluster model, 
hyperspherical harmonic method, lifetime, decay energy. 
PACS: 27.m.+t, 23.50.+Z, 21.10.Tg, 21.10.Dr, 21.60.Gx, 21.45.-v, 21.60.Jz, 
21.10.Jx, 25.70.Mn 



1 Introduction 

In the recent years there were significant advances in the studies of proton 
dripline nuclei. The phenomenon of true two-proton (2p) radioactivity (true 
three-body decay) proposed in Ref. [1] is an important object of these studies. 
The 2p radioactivity is an exclusive quantum-mechanical phenomenon occur- 
ring under specific decay energy conditions {Q2p > 0, Qp < 0) which make 
the sequential proton emission impossible, see Fig. 1. Both protons in this 
case should be emitted simultaneously. The necessary energy conditions are 
widespread near the proton dripline where they are connected to the pairing 
correlations (the effect takes place in isotopes with even mimbcrs of protons). 
A consistent quantum mechanical theory of the 2p radioactivity based on a 
three-body cluster model was for the first time proposed in [2], and further 
developed in [3,4,5,6,7, and Refs. therein]. Being the latest discovered mode 
of radioactive decay (discovered in 2002 in ^^Fe [8,9]) the 2p radioactivity was 
soon found in several other nuclei ^''Zn [10], ^^Ni [11], ^^Mg [13]. The two- 
proton ground-state decays were recently investigated in detail for relatively 
short-lived nuclei ^Be [6] and ^^Ne [13]. The modern trend of this research lies 
in the studies of the proton- proton (p-p) correlations [13,14,6]. These correla- 
tions are expected to be closely connected to the structure of the 2p emitters 
[4,6]. Searches for other 2p-radioactivity cases were proposed in ^''Ar, ^^Ca [15], 
^^Se, ^^'^^Ge, 66,67j^j, i^^j Several experiments dedicated to this phenomenon in 
^Be, ^®Ni, and ^°Ar are now under the way. 

The system ^^S remained somehow "in shadow" while the interest of re- 
searchers were focused on other nuclei. Very little can be found about ^^S 
in literature. The Brookhaven database [16] assigns a half-life of T1/2 ~ 10 
ms and a possible 2p decay mode to this nucleus. Unfortunately, the source 
of this information is not clear. The standard compilation [17] indicates that 
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Table 1 

Energies and charge radii of ^^Si, ^^P, and ^^S nuclei relevant to the 2p-decay prob- 
lem. The energies are given in MeV and the radii are given in fm. The phenomeno- 
logical pairing gap is defined as Apg = 2Qp(^^P) — Q2p{^^S). The corrected 2p-decay 
energy Q2p"(^^S) ~ ^cde — 'S'2n(^^Ne) is based on the calculated Coulomb displace- 
ment energy Ade and on the experimental two-neutron separation energy in the 
mirror system ^^Ne. Maximal and minimal predictions are listed for papers pro- 
viding several sets of results. The two three-body model (3b) results with different 
values are obtained with different charge distributions in 24Si. The results of 
systematics studies are collected in the upper part of the Table and are separated 
from calculations by a horizontal line. 
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Fig. 1. Energy levels and decay scheme for a possible true 2p decay of ^^S. The 
threshold energies are taken from systematics [17]. 

^^S is likely to be nuclear unstable [Q2pC^S) > 0] and is likely to be a true 
two-proton emitter [QpC^^P) > (52p(^^S)], see Table 1. 

The ^^S nucleus has never been studied theoretically in details. It has been 
briefly mentioned in papers [18,19] (phenomenology) as a possible 2p emitter. 
RMF calculations were performed in [20,21,22,23]. The results of these studies 
for are quite controversial (information relevant to and its subsystems 
is summarized in Table 1). The predicted 2p-decay energies are scattered in 
a broad range between and 2 MeV, and the radial characteristics vary con- 
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siderably as well. The overall situation is such that the existing theoretical 
information is not very helpful for planning experimental search for ^®S. It 
is not even clear if ^^S is a true two-proton emitter [Qp{^^P) > Q2p(^^S)] or 
sequential one-proton emission is possible for this nucleus [S'p(^^S) < 0]. The 
lifetime vs decay energy systematics are drastically different in these cases. 
There are theoretical works supporting both possibilities and absolutely no 
experimental data. 

We would hke to provide some support for the true 2p character of the ^^S 
decay using the recent experimental data on ^^Mg [13]. A negative value of 
Sp — —0.42 MeV was inferred for ^^S in Ref. [18] basing on the Ap^ vs Sp 
systematics. Here Apn{Z) = Sn{z^) ~ ^p^n^) is the difference of the nucleon 
separation energies in isobaric mirror partners. It was assumed in Ref. [18] 
that the Apn behaviour should be smooth and thus a value of Ap„ ~ 6 MeV 
was estimated for ^^S (in contrast with Apn — 5.39 predicted by [17]). This 
systematics is shown in Fig. 2 but together with the same ones for the neon 
and magnesium chains. One can easily see that for the known true 2p emitters 
^^Ne and ^^Mg, also belonging to the s-d shell, there is a break in the smooth 
trend of Apn. Such a break could be connected, for example, to structural 
peculiarities of the three-cluster systems in proximity to the proton dripline 
[24]. The Apn value is shifted downwards by 0.35-0.5 McV for the known 
neighbouring systems when we pass the proton dripline. This observation is in 
a good agreement with the prediction of Ref. [17] for supporting the idea 
about the true 2p character of the ^^S decay. 

The discussed situation motivated us to make first a dedicated theoretical 
study of ^^S focusing on the possibility of true two-proton emission and then to 
perform an experimental search for if the theoretical results look promising. 



2 The many-body structure aspect of the 2p decay 



Our approach to the 2p radioactivity problem [2] is based on the approxi- 
mate boundary conditions for the Coulomb three-body problem formulated in 
Ref. [3]. The three-body wave function (WF) ^^^^(p.r^s) with pure outgoing 
boundary conditions and complex energy E = Q2p — iT /2 is constructed for 
an arbitrary small value of F using the hyperspherical harmonic method. In 
this method the WF depends on the collective variable hypcrradius p and on 
the 5-dimensional hyperangle O5. After that the actual F is found using the 
"natural" definition 



F= 4r 



pmax 



N ||^(+)| 



4 



6.0 

^ 5.5 
> 

S 5.0 



4.5 



4.0 



3.5 





m 




;/... 


— • — 


^# 


® 26ri 


------- 

.... ■ .... 


^'Mg - ''Mg 

1 fl-W T 1 Qt. T 




, ♦ ^ 




- "Mg 




-♦ 


1,1, 













3 4 5 
(MeV) 



Fig. 2. Systematics Ap„ vs Sy. For the neon and magnesium chains the values are 
experimental; for the sulphur chain they are partly taken from systematics Ref. [17] 
(dots with circles around). Below the gray dashed horizontal line the value for ^^S 
is positive (it is a true two-proton emitter) and above it is negative (the sequential 
one-proton decay becomes possible). 

The differentials of the flux j defined on the hypersphere with large hyperra- 
dius pmax ~ 10'^ fm and the classical trajectory extrapolation to pext ~ 10^ fm 
are used to define the momentum distributions (see Refs. [4,5,7] for technical 
details). 

The three-cluster approach to the 2p decays is well justified for the closed-shell 
systems or systems with a closed-shell core. However, most of the two-proton 
emitters, including ^^S, do not belong to this class, and important effects of 
the many-body structure can be expected for them. The major components 
of the three-body cluster WFs with = 0^ can be written in a schematic 
spectroscopic notation as 



(+) 



(2) 



The schematic notation [/"]j denotes the Slater determinant of n nucleons 
occupying /j orbital projected on the total spin J and normalized: the three- 
body WF can be considered as normalized in the internal region 



1^ 



(+)||2 



(+)|2 



(3) 



without loss of generality. In papers [4,5] the results of the 2p width calcula- 
tions were provided as functions of weights of the main cluster configurations 
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w{lf) — Xf. Therefore we can write 

nxl)=3{Xl). (4) 

The values which can be put in correspondence with the components of the 
three-body WF (2), to take into account the many-body structure, are over- 
laps of the many-body WFs of the precursor-daughter pair multiplied by a 
combinatorial term 




The overlaps of Eq. (5) are in general case not normalized to unity: 

^2p = EJ^^lVi- (6) 

— 1/2 

If we normalize them to unity, '^a,a-2 '^A,A-2N2p , then the initial many- 
body WF becomes not normahzed to unity ||^a|| = ^2p^- According to Eq. 
(1) the width should be then renormalized as 

r^N,,r{xf/N,,). (7) 

The renormalized width contains a product with structural factor (A'2p) which 
can be easily interpreted as a kind of a spectroscopic factor in analogy with 
the two-body decays. However, contrary to the two-body case, the dependence 
of the width on the structure includes the dependence on the amplitudes Xj. 



They should be adjusted as Xj — )■ Xi/^N2p. That makes the above renor- 
malization a complicated nonlinear procedure. We will see below that the 
many-body effects lead to a significant renormalization of the width beyond 
the Hartree approximation. 

To calculate the amplitudes Xi in a many-body approach we can use a formal- 
ism analogous to the same intrinsic structure model as for direct two-nucleon 
transfer reactions [for instance, pick-up processes {p,t), (n,^He) etc]. In the 
present work we consider the following situation: (i) both the parent nucleus 
Z and the daughter nucleus Z — 2 are essentially many-body systems, (ii) both 
nuclei are spherical in their ground states, (iii) the emitted proton pair trans- 
fers zero total angular momentum, (iv) both nuclei remain in their ground 
states. 



The structure part of the two-nucleon transfer cross section is determined by 
the spectroscopic amplitude B which is defined in the second quantization 
formahsm as [25,26]: 

B(KiKf,kk') = m CjlMiJM{'^Kf\0'li^,)\'i/Ki), (8) 
MMi 
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Fig. 3. (Colour online) Weights of the major components of the overlaps Eq. (5) (left 
axis) and the total overlap normalization as a function of the parameter a (right 
axis). This parameter scales the pairing gap A aA. 



where Ki = {{Ki),Mi} and Kj = {{Kj), Mf} denote quantum number sets 
of the initial and final nuclei with (^Kf, Kf) and without ({Ki)] (Kf)) anerular 
momentum projections Mj, Mj. Index k denotes the full set of the single- 
particle quantum numbers k = {(/c),mfc} in a spherical nucleus. O is the 
two-nucleon transfer operator: 



^i^) ~ 1 + ^ '^jZkh'm,,(^kak' , (9) 

-1- -t- O(kk') rukmy 

where ak is the nucleon annihilation operator in the mean field basis. 

For the ground states of both parent and daughter nuclei we apply the Bardeen- 
Cooper-Schrieffer (BCS) approximation. In this model the ground state WF 
of an even-even spherical nucleus reads: 

\^k) = \K; BCS) = n + ^w44) I-) ' (10) 

k> 

where |— ) is the bare vacuum. In the case of time-reversal symmetry for each 
state > 0, there exists a corresponding time- reversed state < so that all 
the states {fc, k} together form the whole single-particle space. The spectro- 
scopic amplitude defined in Eq. (8) takes the form: 



x^k) = i??,,,^,,,..) = J (jk + - ■ (11) 
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Q2, (MeV) 

Fig. 4. (Colour online) Width of the ■^^S g.s. as a function of the 2p decay energy. 
Solid and dashed curves correspond to the quasiclassical simultaneous emission 
model and to the three-body model, respectively. The RMF-assisted three- 
body model, based on the relativistic Hartree (RH) and on the complete RH+BCS 
results are shown by the dotted and by the dash-dotted curves, respectively. The 
lower hatched area shows the energy and the half-life ranges excluded by the current 
experiment. 

3 Structure of in RMF 

In the present work we employ a RMF model based on the covariant NL3 
density functional [27]. Numerous calculations with NL3 parametrization have 
demonstrated excellent agreement with data in nuclear masses, deformation 
properties [27] and properties of giant resonances and nuclear low-lying re- 
sponse [28]. In RMF model the single- nucleon wave function is a Dirac spinor 
characterized by the set k = {{k),mk} and (k) = {n^, j^, vTfc, ta,} with the ra- 
dial quantum number n^, angular momentum quantum numbers jk, rrik, parity 
TTjt and isospin r^, see Ref. [29] for details. For the pairing interaction we use a 
monopole-monopole ansatz. The BCS gap equation is solved self-consistently 
with the RMF set of equations as described in Ref. [28]. The strength of the 
pairing interaction was adjusted in such a way that the binding energy of ^^Si 
and the phenomenological proton pairing gap in ^^S are well reproduced. 

The results of the RMF calculations important for 2p decay are summarized in 
Table 1 and Fig. 3. The figure demonstrates evolution of the ^^S structure with 
the scaling of the pairing gap. In the relativistic Hartree (RH) case (a = 0) only 
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Fig. 5. (Colour online) Black solid line shows the scaled mean field potential re- 
producing the proton separation energy (right axis). Black dashed line shows the 
potential with a repulsive core excluding the s-wave forbidden state. Single-particle 
s-wave WFs (found on resonance) are shown for these potentials by the blue curves 
of the corresponding style (left axis). 

the nearest to the Fermi surface si/2 single-particle configuration contributes 
to 2p decay and the renormalization of the width is trivial (A^2p = !)• For 
realistic values of the pairing gap (a ~ 1) the valence configuration structure 
for ^^S is characterized by the strong configuration mixing with a dominance 
of the [df\Q configurations. 



4 Three-body 2p decay calculations 

The simple half-life estimates for ^^S are provided in Fig. 4 in terms of the 
models based on systematic values of parameters which we have broadly used 
before. The quasiclassical simultaneous emission model [5] neglecting the inter- 
action between the protons provides a broad "belt" of half-life values between 
the boundaries connected to emission from the [s^]o (upper width limit) and 
[d^jo (lower width limit) configurations. The three-body model [4] predicts 
somewhat narrower belt of possible half-lives (however, it is still as broad as 
two orders of magnitude). 

The RMF results are utilized to modify the three-body calculations in two 
ways, (i) The RMF potentials are used to create the core-p interactions for 
the three-body model, (ii) The three-body internal structure of the ^^S is 
adjusted to be reasonably close to the structure predicted by RMF with a 
help of phenomenological short range potentials in p variable acting on the 
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selected hyperspherical components of the WFs ^3 ' . 

The RMF potential is fitted by the two Woods-Saxon components with Vq — 
{26,-90} MeV, ro = {0.8,2.87}, a = {0.55,0.75} fm (VJ, = -14.3 MeV, 
tq = 2.87, a = 0.75 fm), see Fig. 5. The depth of this potential is scaled to get 
the RMF positions of the p- and (i-wavcs in the nonrelativistic Schrodingcr 
equation. For the s-waves the modified potential is used Vq = {32.7, —19} 
MeV, To = {1.5,4.15}, a — {0.8,0.63} fm, which contains the core and thus 
excludes the Pauli forbidden deep s-state. The 2p width is very sensitive to the 
WF behaviour on the nuclear surface. It can be seen in Fig. 5 that the surface 
behaviour of the WF in the modified potential is practically identical to that 
of the initial potential. The Coulomb potential in the core-p channel is derived 
using the charge distribution of ^^Si obtained in the RMF calculations. The 
simple single Gaussian s-wave potential V{r) — — 31 exp[— (r/1.8)^] is used in 
the p-p channel. 

Thus, the difference of the RMF-assisted three-body results shown in Fig. 
4 compared to results of our ordinary calculations has the following 
three sources, (i) The core-p interactions arc not taken from systematics 
[ro = 1.2(/lc -|- 1)^/^ fm] but are constructed using the RMF results. The 
width is somewhat increased by this choice of the core-p interactions, (ii) The 
three-body structure of ^^S is adjusted to the relative weights of the two-proton 
overlap functions obtained in the RMF calculations. This leads to the width 
reduction compared to the pure three-body calculations, (iii) The 2p width 
value should be renormalized according to the value A^2p of the two-proton 
overlap functions obtained in the RMF calculations. For the predicted struc- 
ture of ^^S this leads to the increase of the width by a factor of two. The latter 
two trends partly compensate each other. Thus, at least for the case of ^^S 
decay, taking the many-body structure of the 2p emitter into account leads to 
the results which are less sensitive to model assumptions. The predictive abil- 
ity of the model is therefore expected to be higher even for quite a simplistic 
method of the structure treatment used in this work. 

This model has important features which we would like to emphasize: (i) the 
internal structure of the three-body WF corresponds well to the structure of 

the overlaps predicted in RMF (ii) penetration through the Coulomb barrier 
is entirely defined by the p-p interaction and the single-particle structure of 
the core-|-p subsystems. 

Predicted complete momentum correlations for the 2p decay of are shown 
in Fig. 6 in terms of the energy distribution parameter e and angle 6k between 
the Jacobi momenta k^,, \s.y: 
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Fig. 6. (Colour online) Kinematically complete correlation patterns for the 2p decay 
of ^^S g.s. in the "T" Jacobi system (e scales the relative energy between two 
protons) and in the "Y" Jacobi system (e scales the energy between the daughter 
nucleus ^^Si and one of the protons). Angles 9k are defined in Eq. (12). 



e = E^/Et , cos{9k) = {K ■ ky)/(fcx ky) , 
ET = E, + Ey = kl/2M^ + kl/2My , 
^ Aaki - Aika ^ A3(ki + ka) - {Ai + ^2)1^3 



where and My are the reduced masses of the X and Y subsystems (see, 
e.g. [6] for details). The correlation patterns demonstrate the regions of strong 
suppression, connected with the repulsive core-proton Coulomb interaction 
[e ~ 0.5, cos(6'fc) ~ ±1 in the "T" system]. Another feature is a broad peak 
[cos(^fc) ~ in the "T" system] connected with the large weight of the s-wave 
configuration in ^^S. Similar correlation patterns have been obtained for other 
s/d-she\\ 2p emitters [4] and were later found to be in agreement with the 
experimental data in Ref. [31]. 

Results of the theoretical studies indicate that the considerations for further 
experimental searches for are especially complicated. If we take the stan- 
dard energy prescription made by Audi and Wapstra [17] then we have a 50% 
"probability" that can be studied in an implantation experiment (half-life 
exceeding ~ 100 ns - the typical flight-time for the fragment separation). Oth- 
erwise the decay-in-flight technique [13] can be used for the ^^S study down 
to half-life times as short as some picoseconds. The Q2p predictions made by 
different theoretical models point at higher energies Q2p ^ 0.8 MeV but do 
not seem to provide a solid basis for the choice of experimental technique. We 
decided to make the first attempt to study ^^S in the implantation experiment. 
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Fig. 7. (Colour online) Ai^-ToF identification plots for the cases of separator tuning 
on (a) and (b). 



5 Experiment 



The experiment was performed at the ACCULINNA fragment separator [30] 
(JINR, Dubna, Russia). A 50.3 A MeV beam of ^2311+ ions delivered by the 
U400M cyclotron with intensity up to 30 pnA bombarded a 92.4-mg/cm2 
beryllium production target installed in the separator object plane. A 92.4- 
mg/ cm2 beryllium wedge was installed in the intermediate focal plane of the 
separator. Combined with a slit restricting the momentum acceptance to about 
2.0%, this rather thick wedge reduced in number the set of ^23 fragmentation 
products arriving at the achromatic focal plane when the separator was tuned 
for a desired product with its specific charge Z and mass A. Mostly, the nuclei 
arriving in the beam cocktail at the final focal plane of the separator were 
severely depleted in intensity. This provided favourable conditions to detect 
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^^S anticipated to be a very rare product of the ^^S fragmentation. It turned 
out that only a few fragmentation products with different Z and A occurred 
to be close in transmission probability to the searched nucleus. This appeared 
to be useful as such "satellites" were used for tuning the separator when rare 
fragments, such as ^^S, were searched for. 

To ensure unequivocal identification for the nuclei passing through the sepa- 
rator we used both the time-of-flight-A£^ (ToF-A£') and IS.E-E techniques. 
The ToF array, having about 0.6 ns time resolution, was placed on a 8.5-meter 
flight base between the achromatic and final focal planes of the separator. The 
ToF detectors were two 125-/im thick BC-418 plastics inclined at 45° to the 
separator axis and each plastic viewed by its pair of XP2020 photomultipliers. 
The telescope composed of one 68-/xm and one 1000-/xm, 6x6 cm^ silicon strip 
detectors (8x8 strips) was installed in the final focal plane. Besides the AE-E 
identification the hit positions of the nuclei were measured by this telescope. 
Also, the signals from the telescope AE detectors were used for the ToF-A£^ 
identification. In all the runs the load of the strip detectors did not exceed 
500 pps. Taking into account the granularity of these detectors such counting 
rate makes pileups negligible. 

Four runs were performed in which the separator was sequentially tuned to 
give the optimum transmission for ^^S, ^^S, ^''S and ^^S. The measurements 
took 870, 1200, 4500 and 1.1 x 10^ seconds, respectively. For the two latter 
cases the AE'-ToF plots are presented in Fig. 7. No events of ^^S, as well as 
^^P, were found in the measurements. To extract the upper lifetime limits for 
these nuclei, the isotope yields of some known nuclei from sulphur to neon 
were systematized depending on their reaction Q values, see Fig. 8. In this 
figure the results of runs tuned for the sulphur isotopes are shown. In addition 
to the three sulphur isotopes (^^S, ^^S and ^^S) the yields of nuclei represent- 
ing the "sateUites" in the four separator tunes are presented in Fig. 8 (a). 
These satellites are the nuclei having parameter A^^^/Z^^^ close to that of 
the tuned sulphur isotopes (variation from —4% to +3% was allowed). This 
parameter defines the selectivity of the energy loss achromat. For the nuclei 
selected in such a way the yield dependence on the reaction Q values should 
have exponential character reasonably close to the exponential dependence 
of the production cross sections. Both the yield and the cross section depen- 
dence on ^reaction IS givcu in Fig. 8 (b) for a subset of these nuclei having 
high transmission rates through the fragment separator, namely above 0.25%. 
Naturally, these appear to be sulphur, phosphorus, and silicon isotopes. The 
cross sections were defined using the transmission rates calculated with the 
use of relations from Ref . [32] . 

The main contribution to the experimental errors stems from the systematic 
uncertainties of both the fiux of the ^^S primary beam and the transmission 
through the separator. It was estimated to be about 25%. The extrapolation 
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of the obtained systematics shows that in the presented experiment one could 
expect the observation of I6I4 ^iid SGOtyo events of and ^^P, respectively. 
The total length of the ACCULINNA separator from the production target 
to the final focus position is 21.3 m. The full times of flight for these nuclei 
from the production target to the detector system were 314 ns and 322 ns 
(accounting for energy losses in the wedge and ToF-plastics). Using the Poisson 
distribution the half-life limit for can be found as Ti/2(^^S) < 79 ns with 
63% confidence level and Ti/2(^^S) < 157 ns with 98% confidence level. In 
these estimations one and four events respectively were assumed instead of 
observed zero events. The 79 ns half-life limit is in qualitative agreement with 
nonobservation of ^^S as an expected accompanying nucleus in the recent 
experiment with a ^^Si secondary beam at MSU [33]. 

The corresponding half-hfe hmits for ^^P are found as Ti/2(^^P) < 38 ns and 
7i/2(^^P) < 50 ns. The 38 ns half-life limit is compatible with the limit of 30 
ns provided in Refs. [16,17,34] (this value seems to be based on an unpublished 
work). It can be easily estimated that such a limit provides quite a relaxed 
limit for the proton decay energy in ^^P: Qp{^^P) ^110 keV (according to 
an R-matrix estimate with systematic parameters). This is not a restrictive 
result to define the character of the ^^S decay (true or sequential 2p decay, 
see Table 1 and relevant discussion in the introduction). However, it helps to 
provide a two- proton decay energy limit for ^^S. 

Using the dependence shown in Fig. 4 we find that a half-hfe hmit of 79 ns 

provides a 2p-decay energy limit <52p(^^S) > 640 keV assuming the true two- 
proton decay. If the one-proton decay were possible for this nucleus (which we 
find to be unlikely situation, but cannot exclude completely) one would get 
a hmit for one-proton decay energy Qp(^'^S) > 120 keV (R-matrix estimate). 
Combining it with Qp(^^P) hmit we obtain a hmit (52p(^^S) > 230 keV. The 
previously reported half-life value of 10 ms [16] belongs to the timescale where 
the decay can be attributed to the weak transitions as well (nucleus with such 
a half- life could appear to be even nuclear stable). The half-life times in the 
nanosecond range and shorter would undoubtfuUy identify the ^^S nucleus as 
a 2p emitter. 



6 Conclusions 

We have studied both theoretically and experimentally properties of the pro- 
ton dripline nucleus ^^S. Structures of ^^S and its subsystems (^^Si, ^^P) were 
obtained in the RMF approach. The 2p decay properties of ^^S were obtained in 
a three-body model tuned by the RMF results. It is demonstrated that in this 
joint approach the results appear to be noticeably different than the results of 
the three-body "P" model based on a systematics input. The uncertainties of 
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Fig. 8. (Colour online) Panel (a) shows the yields of different nuclei as a function of 
the reaction Q value. In this panel the symbol size is comparable with the error bar 
values. Panel (b) shows the yields (lower group of symbols, red line) and production 
cross sections (upper group of symbols, blue line) for the isotopes with highest 
transmission rate through the fragment separator. 

the 2p width connected with structure of 2p emitters are significantly reduced 
in this approach. 

Production of ^^S in the fragmentation reaction was studied experimentally. 
We have found that either the production of ^^S is anomalously small (more 
than one order of magnitude than expected according to systematics) or the 
lifetime is too short. The systematical value of the production cross section 
enables us to estimate the upper half-life limit of ^^S as 79 ns. Theoretical 
lifetime dependence on the 2p-decay energy allows one to infer the lower limit 
of 640 keV for the Q2p of ^^S. This value is in principle consistent with the 
theoretically predicted range of mass excess for ^^S but it narrows significantly 
its uncertainty. In the case that ^^S has an open one proton decay branch (se- 
quential two-proton decay) a less restrictive limit Q2p > 230 keV is obtained. 

Our results are in contradiction with the existing information about ^^S as a 
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possible 2p emitter with the half-hfe of about 10 ms [16]. Our results allow 
focusing the search for ^^S in the half-life range 10 ps - 100 ns where the 
decay-in-flight technique [13] is the adequate experimental method. 
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